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The  radiation  fog  model  ia  being  developed  in  an  attempt 
to  predict  the  initiation,  development  and  characteristics 
of  these  fogs.  The  model  is  primarily  baaed  upon  three 
Important  processes  which  interact  to  determine  whether  a 
fog  will  form  in  preference  to  haze  development  and  what 
characteristics  and  fluctuations  are  likely  to  occur  if  a 
fog  does  form.  These  processes  are: 

1.  Radiative  exchanges' 

2.  Droplet  growth  and  CCN  spectra’ 

3.  Turbulence  characteristics. 

The  exact  nature  of  these  interactions  is  poorly  under¬ 
stood  and  consequently  have  not  yet  been  fully  included  in 
existing  models.  To  date,  therefore,  accurate  simulation  of 
all  fog  characteristics  on  various  temporal  and  spatial  scales 
has  not  been  achieved. 

Data  from  recent  field  studies,  especially  Meppen,  have 
produced  results  which  will  enable  not  only  the  model  to  be 
fully  tested,  but  also  greatly  Bssist  in  the  understanding  of 
radiation  fog,  which  will  in  turn  allow  more  accurate  and 
detailed  models  to  be  developed. 
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MODEL  FORMULATION 

The  model,  at  present,  utilises  the  following  equations 
to  describe  the  changes  in  temperature,  supersa turation ,  and 
liquid  water  content,  in  addition  to  explicitely  calculating 
droplet  growth. 


Droplet  growth 
fir. 


6  t 


“■Wki 


B 


CM. 


2  +  l  -  DFQ  (r) 

U 


-] 


where  the  coefficients  are 
A 


lhPl  ^lhm  PlV 

=  kt  V«gt  ■  J  + 


m 


B 


=  2p  M 

Plrgt 


C  =  6.92L6  X  10"5m3kg-1 


°  =  -0 


Temperature 

1  £  ♦  fr  (-H )  ♦  ijfa 


il 

6  t 


P bCP  fiZ 


Super saturation 
1 


jgS 

fit 


Mr 


6MR  _  LhM 


RS  6 1 
Liquid  Deter  Content 


_  6  T 

RqT2  6 1 


AM  _  i 


£ —  Ij  Z!  i\j  r 


fit 


3  PD  fit  i  j  "ijrij 


Condensation  rate 


cn  -  St  -  fl  -  h  (1i) 


Mixing  ratio 


Me 


(1  +  S)  M 


RS 


(1) 


(2) 


(3) 


(A) 


(5) 


V  V  *,*  V  *. 

i.fr  alkmTa  .fa  V 


T  x  dry  bulb  temperature 
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pf  =  density  of  wBter 
k  =  thermal  conductivity  of  Bir 
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s  molecular  weight  of  water 
x  universal  gaa  constant 
x  diffusivity  of  water  vapour  in  air 
x  surface  tension  of  water  vapour 
=  density  of  air 
x  net  radiative  flux 
=  turbulent  exchange  coefficient 
x  potential  temperature 
x  latent  head  of  vaporization 
x  condensation  rate 

x  specific  heat  of  air  at  constant  pressure 
x  saturation  mixing  ratio 


RS 

MR  =  mixing  ratio 

r^  =  droplet  radius  class  i 

nij  x  nucleus  mass 

x  concentration  of  drops  in  class  i 
F  x  radiative  exchange  function  for  droplet  growth 

SBC 
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Q  (r)  x  absorption  efficiency  Df  water  droplets 
8 

x  1.18(1  -  exp(-0.28r^ j)> 

111  x  liquid  water  content 

G  =  gravitational  sedimentation 

SBC  x  Stef an-Bol tzmann  constant 

y  =  function  of  terminal  velocity  of  drops 


The  droplet  growth  is  determined  by  Eq.1,  which  in 
addition  to  the  usual  terms,  contains  a  term  derived  by 
Roach  (1976)  which  allows  for  net  radiative  loss  from  the 
droplet.  Eq.  2  describes  the  time  rate  of  change  of  temperature 
which  is  determined  by  the  divergence  of  net  radiation,  the 
eddy  diffusion  of  heat,  and  the  release  of  heat  through  con¬ 
densation.  The  change  of  supersaturation ,  eq.  3,  is  a  result 
of  changes  in  temperature  and  removal  by  condensation. 


In  addition  to  these  equations,  changes  in  droplet 

— 

concentration,  determined  by  eddy  diffusion  and  gravitational 


sedimentation,  are  also  included. 
Number  concentration 


Y 


Gravitational  sedimentation 


Although  these  terms  are  included  in  the  model  their 
formulation  (and  to  a  large  extent  that  of  other  terms) 
depend  on  the  trubulence  characteristics  of  fogs  which  are 
not  well  known.  It  is  believed  that  the  value  of  K,  the 
turbulent  exchange  coefficient,  and  its  method  Df  derivation, 
are  critically  important  in  determining  the  precise  nature  of 
fluctuations  in  the  fog,  interactions  between  the  properties 
of  the  fog  top  and  the  fog  structure  below,  and  interactions 
between  the  fog  and  the  ground.  The  exact  nature  of  these 
interactions  are  being  closely  examined. 
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INPUT  PARAMETERS 

The  Input  parameters  for  the  model  are  88  follows: 

1.  Mixing  ratio  -  assumed  constant  with  height.  In  addition 
to  temperature  determines  initial  relative  humidity. 

2.  Isothermal  atmosphere  up  to  AOOm  -  the  ground  temperature 
is  specified. 

3.  Turbulent  exchange  coefficient  -  this  can  be  varied  as  an 
input  parameter  but  also  can  be  re-calculated  according 
to  variations  in  stability. 

U.  Soil  temperature. 

5.  Net  radiative  flux  -  appropriate  values  can  be  input  or 
calculated  from  a  parameterization  scheme. 

6.  CCN  spectra  -  the  number  concentration  Df  nucleus  sizes 

ere  input,  but  the  model  assumes  them  to  be  soluble  ammonium 
sulphate  ((NH^JgSO^).  However,  the  chemical  composition 

of  the  CCN  can  be  varied  as  required.  The  aerosol  concentration 

is  in  the  form  of  a  Junge  distribution 
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MODEL  PROCEDURE 

-  i 

Droplet  growth  (b  brief  synopsis) 

At  first,  while  relstive  humidity  is  low,  the  model  is 
run  with  a  U  sec  time  step,  although  this  can  be  changed  as 
necessary.  A  shorter  time  step  would  be  desirable,  but  com¬ 
putation  time  would  be  much  greater  with  little  improvement 
in  accuracy  in  the  present  model. 

At  first  droplets  are  set  at  their  equilibrium  radii, 
which  are  determined  by  nucleus  mass,  chemical  composition, 
and  relative  humidity,  and  are  calculated  in  a  subroutine 
of  the  program.  The  initial  stage  of  the  model  can  be  used 
therefore  to  describe  reasonably  well  the  development  of 
haze  droplets  before  fog  initiation. 

At  each  time  step  in  the  model  the  equilibrium  radii 
ere  calculated,  so  as  allowing  an  upper  limit  to  growth  to  be 
initially  specified.  Small  droplets,  with  a  U  sec  time  step, 
reach  equilibrium  radius  within  that  time  step,  but  as  relative 
humidity  and  drop  size  increase,  so  the  larger  drops  fail  to 
attain  their  equilibrium  size. 

When  the  relative  humidity  in  the  lower  layers  has 
exceeded  99.99%,  then  the  time  step  is  reduced  to  1  sec  for 
more  sensitive  computation.  The  droplets  continue  to  grow 
in  response  to  the  increase  in  relative  humidity  and  finally 
supersaturation.  The  radiative  loss  term  in  the  droplet  growth 
equation  allows  some  drops  to  be  activated  before  100%  relative 
humidity,  by  lowering  their  critical  supersaturation.  Essentially 
this  radiation  term  describes  the  reduction  of  droplet  temper¬ 
ature  in  relation  to  environmental  temperature,  and  allows 


some  drops  to  continue  to  grow  in  an  undersaturated  atmosphere. 

The  procedure  for  calculating  the  other  parameters  required 
for  fog  development  will  be  set  out  fully  in  the  program  docu- 
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mentatiun  86  explanation  in  the  context  of  this  report  is 
difficult  without  referring  extensively  to  a  procram  descrip¬ 
tion.  A  more  detailed  description  of  droplet  growth  will  also 
be  provided  with  the  program  documentation. 
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